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^ . Abstract 

A distributed detection scheme where the sensors transmit with constant modulus signals over a 
Gaussian multiple access channel is considered. The deflection coefficient of the proposed scheme is 
shown to depend on the characteristic function of the sensing noise and the error exponent for the 
system is derived using large deviation theory. Optimization of the deflection coefficient and error 
exponent are considered with respect to a transmission phase parameter for a variety of sensing noise 
distributions including impulsive ones. The proposed scheme is also favorably compared with existing 
amplify-and-forward and detect-and-forward schemes. The effect of fading is shown to be detrimental 
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to the detection performance through a reduction in the deflection coefficient depending on the fading 
statistics. Simulations corroborate that the deflection coefficient and error exponent can be effectively 
used to optimize the error probability for a wide variety of sensing noise distributions. 

Index Terms 

Distributed Detection, Multiple Access Channel, Constant Modulus, Deflection Coefficient, and 
Error Exponent. 

I. INTRODUCTION 

In inference-based wireless sensor networks, low-power sensors with limited battery and peak 
power capabilities transmit their observations to a fusion center (FC) for detection of events or 
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estimation of parameters. For distributed detection, much of the literature has focused on the 
parallel topology where each sensor uses a dedicated channel to transmit to a fusion center. 
Multiple access channels offer bandwidth efficiency since the sensors transmit over the same 
time/frequency slot. The sensors may perform local detection through quantization in which 
case the decision is encoded into a specific waveform to be sent to the FC. Instead, sensed 
information may be sent to the FC using analog modulation which transmits the unquantized data 
by appropriately pulse shaping and amplitude or phase modulating to consume finite bandwidth. 

In [1], the distributed detection over a multiple access channel is studied where arbitrary 
number of quantization levels at the local sensors are allowed, and transmission from the sensors 
to the fusion center is subject to both noise and inter-channel interference. References [2]- 
[5] discuss distributed detection over Gaussian multiple access channels. In [2], detection of a 
deterministic signal in correlated Gaussian noise and detection of a first-order autoregressive 
signal in independent Gaussian noise are studied using an amplify-and-forward scheme where 
the performance of different fusion rules is analyzed. In [3], a type-based multiple access scheme 
is considered in which the local mapping rule encodes a waveform according to the type [6, pp. 
347] of the sensor observation and its performance under both the per-sensor and total power 
constraints is investigated. This scheme is extended to the case of fading between the sensors 
and the FC in [4] and its performance is analyzed using large deviation theory. In the presence of 
non-coherent fading over a Gaussian multi-access channel, type-based random access is proposed 
and analyzed in [5]. In [7], the optimal distributed detection scheme in a clustered multi-hop 
sensor network is considered where a large number of distributed sensor nodes quantize their 
observations to make local hard decisions about an event. The optimal decision rule at the 
cluster head is shown to be a threshold test on the weighted sum of the local decisions and its 
performance is analysed. 

Two schemes called modified amplify-and-forward (MAF) and the modified detect-and-forward 
(MDF) are developed in [8] which generalize and outperform the classic amplify-and-forward 
(AF) and detect-and-forward (DF) approaches to distributed detection. It is shown that MAF 
outperforms MDF when the number of sensors is large and the opposite conclusion is true 
when the number of sensors is smaller. For the MDF scheme with identical sensors, the optimal 
decision rule is proved to be a threshold test in [9]. Decision fusion with a non-coherent fading 
Gaussian multiple access channel is considered in [10] where the optimal fusion rule is shown 
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to be a threshold test on the received signal power and on-off keying is proved to be the optimal 
modulation scheme. A distributed detection system where sensors transmit their observations 
over a fading Gaussian multiple-access channel to a FC with multiple antennas using amplify- 
and-forward is studied in [11]. In all these cases, the sensing noise distribution is assumed to be 
Gaussian. Even though the Gaussian assumption is widely used, sensor networks which operate 
in adverse conditions require detectors which are robust to non-Gaussian scenarios. Moreover, 
in the literature there has been little emphasis on distributed schemes with the desirable feature 
of using constant modulus signals with fixed instantaneous power. 

A distributed estimation scheme where the sensor transmissions have constant modulus signals 
is considered in [12]. Distributed estimation in a bandwidth-constrained sensor network with a 
noisy channel is investigated in [13] and distributed estimation of a vector signal in a sensor 
network with power and bandwidth constraints is studied in [14]. The estimator proposed in [12] 
is shown to be strongly consistent for any sensing noise distribution in the iid case. Inspired 
by the robustness of this estimation scheme, in this work, a distributed detection scheme where 
the sensors transmit with constant modulus signals over a Gaussian multiple access channel is 
proposed for a binary hypothesis testing problem. The sensors transmit with constant modulus 
transmissions whose phase is linear with the sensed data. The output-signal-to-noise-ratio, also 
called as the deflection coefficient (DC) of the system, is derived and expressed in terms of the 
characteristic function (CF) of the sensing noise. The optimization of the DC with respect to the 
transmit phase parameter is considered for different distributions on the sensing noise including 
impulsive ones. The error exponent is also derived and shown to depend on the CF of the sensing 
noise. It is shown that both the DC and the error exponent can be used as accurate predictors of 
the phase parameter that minimizes the detection error rate. The proposed detector is favorably 
compared with MAF and the MDF schemes developed in [8], [9] for the Gaussian sensing noise 
and its robustness in the presence of other sensing noise distributions is highlighted. The effect 
of fading between the sensors and the fusion center is shown to be detrimental to the detection 
performance through a reduction in the DC depending on the fading statistics. Different than [12] 
where the asymptotic variance of an estimator is analyzed, the emphasis herein is on derivation, 
analysis, and optimization of detection-theoretic metrics such as the DC and error exponent. Our 
aim in this paper is to develop a distributed detection scheme where the instantaneous transmit 
power is not influenced by possibly unbounded sensor measurement noise. 
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The paper is organized as follows. In Section II, the system model is described with per-sensor 
power constraint and total power constraint. In Section III, the detection problem is described 
and a linear detector is proposed. The probability of error performance of the detector is analyzed 
in Section IV. The DC is defined and its optimization for several cases is studied in Section V. 
The presence of fading between the sensors and the fusion center is discussed in Section VI. 
The error exponent of the proposed detector is analyzed in Section VII. Non-Gaussian channel 
noises are discussed in Section VIII. Simulation results are provided in Section IX which support 
the theoretical results. Finally, the concluding remarks are presented in Section X. 

II. System Model 

Consider a binary hypothesis testing problem with two hypotheses H , Hi where P , Pi are 
their respective prior probabilities. Let the sensed signal at the i th sensor be, 

{9 + rii under Hi 
(1) 
rii under H 

i = 1,. . . ,L, # > 1 is a known parameter whose presence or absence has to be detected, 
L is the total number of sensors in the system, and rii is the noise sample at the i th sensor. 
The sensing noise samples are independent, have zero median and an absolutely continuous 
distribution but they need not be identically distributed or have any finite moments. We consider 
a setting where the i th sensor transmits its measurement using a constant modulus signal ^fpe^ Xi 
over a Gaussian multiple access channel so that the received signal at the FC is given by 

L 

i=i 

where p is the power at each sensor, cu > is a design parameter to be optimized and v ~ 
CJ\f(0, &l) is the additive channel noise. We consider two types of power constraints: Per-sensor 
power constraint and total power constraint. In the former case, each sensor has a fixed power 
p so that the total power P T = pL, and as L — > oo, P T — > oo; in the later case, the total 
power P T is fixed for the entire system and does not depend on L, so that the per-sensor power 
p = P T /L ->■ as L ->■ oo. 

'the proposed scheme will work without any difference for 9 < due to symmetry if we substitute —9 in the place of 9 in 
all the equations. 



November 16, 2010 



DRAFT 



5 



III. The Detection Problem 



The received signal y L under the total power constraint can be written as 

L 



l~P~ L 

y L = J^Y, ej " Xi + v - ^ 




i=i 

We assume throughout that P = Pi = 0.5 for convenience even though other choices can be 
easily incorporated. With the received signal in (3), the FC has to decide which hypothesis is 
true. It is well known that the optimal fusion rule under the Bayesian formulation is given by: 

nmm gg> = 1 (4) 

f(y L \H ) h Pi 

where f(y L \Hi), is the conditional probability density function of y L when Hi is true. The 
equation (3) can be rewritten as follows: 

. ^ —(ujXi) 1 + j d — I sm(uxi) I + v. 

Since there are L terms in the first summation involving the cosine function, we need to do L 
fold convolutions with the PDFs of cos(wrrj) and another set of L fold convolutions with the 
PDFs of sin(cuxj). Then we need to find the joint distribution of the PDFs obtained thus for the 
cosine and sine counterparts. This joint PDF will need to be convolved with the PDF of v. It 
is not possible to obtain a closed form expression for these (2L + 1) fold convolutions. Hence, 
f(y L \Hi) is not tractable. Therefore, we consider the following linear detector which is argued 
next to be optimal for large L: 

U[y L e-^] - n VL } I , (5) 

where we define 3?[y] as the real part, and Q[y] as the imaginary part of y. Note that the detector 
in (5) would be optimal if y L were Gaussian. Clearly due to central limit theorem y L in (3) 
is asymptotically Gaussian, which indicates that (5) approximates (4) for large L. With the 
Gaussian assumption, the variances of y L in (3) under the two hypotheses are the same and 
given by Var(y L \H ) =Var(y L \H 1 ) = [Pt(1 — V^t^)) + a v\> wnere fn{w) is the characteristic 
function of rij. Hence, the optimal likelihood ratio simplifies to the detector in (5) which is linear 
in y L , when y L is assumed Gaussian which holds for large L. However as will be seen in Section 
IV, we do not assume that y L is Gaussian for any fixed L when we analyze the performance 
of the detector in (5) or in finding the associated error exponent in Section VII. We proceed by 
expressing the probability of error. 
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IV. Probability of Error 
The detector in (5) depends on the design parameter to and this means that the probability of 



error will in turn depend on u. Let P e (u) be the probability of error at the FC: 

Po(tu) 



- Pr [error|iJo] H — P r [error] if i] = Pr [error \H ] 
2 2 



(6) 



where Pr [error \H^\ is the error probability when Hi, i e {0,1}, is true and the last equality 
holds due to symmetry between the two hypotheses which is explained as follows. From the 
detection rule (5), the probability of error under H is given by 



Pr [error \H ] = Pr [X\y L ] < X\y L e- ju& \ \H ] , 
where the received signal in (3) under H is given by 

<)l = [ Vj;Yl cos ( wn i) + V jE sin(wni) + 3f[u] 

i=l / V i=l 

Substituting (8) for y L in (7) and doing some algebraic simplifications we get, 



(7) 



(8) 



Pr [error | H ] = Pr 



E 

i=l 

V 



2 sin 



c< is | wrii - — + - 



+ 




v T < 



(9) 



where t> T := 3?[t>](l — cos(cu9)) — $s[v]sm(uj9). Similarly, Pr [error \Hi] is same as that of 
(9) except the argument of the cosine function is replaced by (urn* + uO/2 — n/2). To see the 
symmetry between the two hypotheses asserted in (6), let ( := (cu9/2 — ir/2) for convenience, so 
that cos(wnj=FC) — [cos(cunj) cosC + sin(±cunj) sin^]. Since n» is symmetric, wn^ and — oorii have 
the same distribution which implies that the random variables cos(amj — () and cos(wnj + C) have 
the same distribution establishing that Pr [error |ifi] = Pr [error|if ]- Therefore, the probability 
of error in (6) is given by (9). We are interested in using (9) to find the cu that minimizes 
the probability of error at FC. Since P e (ui) is not straightforward to evaluate, we optimize two 
surrogate metrics to select to. These are the error exponent and the DC. The error exponent is an 
asymptotic measure of how fast the P e (co) decreases as L — > oo, and is specific to the detector 
used in (5) and will be considered in Section VII. The DC, on the other hand, is specific to the 
model in (3), and does not depend on any detector. 
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V. Deflection Coefficient and its Optimization 

We will now define and use the deflection coefficient which reflects the output- signal-to-noise- 
ratio and widely used in optimizing detectors [15]— [18]. The DC is mathematically defined as, 

D( U) ~ i MLz^W . (10) 

L vaT[y L \H \ 

By calculating the expectations in (10), it can be easily verified that the DC for the signal model 
in (2) is given by: 

= 2^M[l-cosM)] (U) 



where (p n (oj) = E[e? UUi ] is the CF of rij. The CF (p n (u>) does not depend on the sensor index i, 
since we will be initially assuming that rii are iid. We will consider the non-identically distributed 
case in Section V-D. Note that D(oo) > and that (p n (u>) is real-valued since is a symmetric 
random variable. Moreover, (p n (oj) = (p n (—u>) so that D(ui) = D(—oo) which justifies why 
we will focus on oo > throughout. The factor (l/L) introduced in (10) does not appear in 
conventional definitions of the DC but included here for simplicity since it does not affect the 
optimal oo. 

A. Optimizing D(uo) 

We are now interested in finding u by optimizing D(u>): 

oo* := argmaxll(w). (12) 

Since f n {oS) < 1, when al > 0, D(u) is bounded, and achieves its smallest value of D(uS) = 
as uj — > 0. On the other hand, as u — > oo, we have lim w ^ 00 D(w) = 0. This implies that the 
maximum in (12) cannot be achieved by oo = or oo = oo and establishes that there must be a 
finite oo* E (0, oo) which attains the maximum in (12). 

In what follows, we will further characterize oo* by assuming that (p n (oo) > and tp' n {oo) < 
for all oo > 0. Many distributions including the Laplace, Gaussian and Cauchy have CFs that 
satisfy this assumption. Indeed all symmetric alpha-stable distributions [19, pp. 20] of which the 
latter two is a special case, satisfy this assumption. We now have the following theorem which 
restricts oo* in (12) to a finite interval. 

Theorem 1. If <p n (u>) is decreasing and differentiate over oo > 0, then oo* e (0, ir/6). 
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Proof: First, note that (p n (u) > which is implied by the assumption that (p n (uj) is 
decreasing and that ip n (uj) — > as u — > oo. Let D(u) = C(u)[l — cos(w$)] with C(u) := 
2(Pn(u})/[l — <Pn(uj) + <?v/Pt] for brevity. Since (p n (uj) is decreasing onw>0 and <p n {u) > 0, 
C(oo) is also decreasing. Because [1 — cos(cu#)] is periodic in u with period 2n/9, 

D ( u + % ) = [1 - cos(u;0)]C ( w + % j < [1 - cos(u;0)]C(u;) = (13) 



/ V , 

Noticing that D(2n/9) = which rules out u* = 2n/9, we have u* e (0,2n/9). To further 
reduce the range of uo* by half, consider the fact that D(0) = D(2tt/9) = 0, which combined 
with D(u) > for u G (0,2n/9) implies that oo* e (O,2?r/0) satisfies £>'(w*) = 0. Writing 
D'(oj*) =0 we obtain: 

[esin(u*e)] gV) 

[cos(w*0)-l] C(w*)' 
Since C(u;) > is decreasing, the right hand side (rhs) of (14) is negative and it follows that 
oj* e (Q,ir/0) as required. ■ 
By the definition of to*, it is clearly a function of 9. We showed in Theorem 1 that < to* < 
n/9 if f' n {oS) < for u; > 0. Note that when u — 0, there is no phase modulation done, and what 
is transmitted is a constant signal which actually contains no information about Xj. Therefore 
the boundary value w = is not a valid choice. When uo = n/9, the detector in (5) actually 

Ho 

simplifies to: 3?[yzJ ^ 0. While u; = 7T/6 1 is a valid choice, it is optimal only when 9 is large as 

Hi 

will be proved in Theorem 2. We now investigate the behavior of u* when 9 is large without 
assuming anything about <p n {u) except the absolute continuity of its distribution, and show that 
uo* ps -n/9 for large 6> in the sense that cu*9 — > 7r, as 9 — > oo. 

Theorem 2. If a 2 v > 0, and are ad an J /zave absolutely continuous distributions, 

lim w*0 = 7T. (15) 

0— »0O 



/ 7T \ 4Pt 

^ ( ^ ) < ^(w*) < sup[l - cos(w0)] sup C(w) = — ^- , (16) 

\9/ uj>Q uj>0 <J„ 



Proof: We have 

cv) = 

where the first inequality is because maximizes D(u), and the second inequality follows from 
-D(u;) = C(u})[l — cos(u9)]. Recalling that lim w ^ (p n (ui) = 1 we take the limit as — >• oo in (11) 
and obtain Hindoo D(ir/0) = AP T /a 2 , which using (16) shows that Hindoo D{u*) = 4P T /a 2 . 
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Since ip n (0) > <p n (u) and because D(u>) is an increasing function of ^(cu), from (11) it is clear 
that the only way Hindoo D{uj*) = AP^/crl holds is if oo* — > and oo*9 — > tt, as 9 — > oo. ■ 
Theorem 2 establishes that when 9 is large we have an approximate closed-form solution for 
oo* ps 7T/6 1 for any absolutely continuous sensing noise distribution. 

5. Finding the Optimum oo for Specific Noise Distributions 

Theorem 1 showed that oo* £ (0, ir/9) for a general class of distributions. Under more general 
conditions, Theorem 2 establishes that oo* ps n/ 9 when 9 is large. To find oo* exactly, we need 
to specify the sensing noise distribution through its CF, (p n (oo). In what follows we describe 
how to find oo* for several specific but widely used sensing noise distributions. We will assume 
throughout that the assumptions of Theorem 1 ((p' n (co) < for oo > 0) are satisfied so that 
oo* £ (0,ir/9), which holds for Gaussian, Cauchy and Laplacian distributions, among others. We 
will assume al > throughout this subsection. 

1) Gaussian Sensing Noise: In this case, we have (p n (oo) = e~ w2 °™/ 2 so that <Pn(oo) = e~ u)2(J ™, 
where is the variance of n^. To simplify (11) we substitute (3 = oo9. Since oo £ (0,7r/6>) we 
have (3 £ (0,7r). Note that the value of oo that maximizes (11) over oo is related to the (3 that 
maximizes D(j3/9) through the relation oo = [3/9. Differentiating D((3/9) with respect to (3, 
equating to and simplifying we obtain, 



with a :— [1 + (<r 2 /P T )]- Equation (17) can not be solved in closed-form. However it does have 
a unique solution in (3 £ (0, ix) as shown in Appendix 1. The solution to (17), /3q, can be found 
numerically and the optimum oo for the Gaussian case is Uq = (3q/9. 

2) Cauchy Sensing Noise: In this case, (p n (oo) = e~ lw so that ^(uo) = e~ 2lw where 7 is the 
scale parameter of the Cauchy distribution. It is well known that no moments of this distribution 
exists. Substituting tp n {oo) in D(oo) and letting (3 = oo9 we have, 



with a := [1 + (<7^/P T )] and (3 £ (0,7r). It can be verified that the equation (18) has a unique 
maximum over (3 £ (0,7r) as shown in Appendix 2. The /?£ that maximizes (18) can be found 
numerically and 00 q = (3q/9. 




(17) 




(18) 
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When c 2 /Pt is sufficiently large (i.e., the low channel SNR regime) compared to [1 —(p 2 n {u))] 
in D(oj), the problem in (11) can be transformed into maximizing (p 2 n (u)[l — cos(ui6)] over 
cj G (0, tt/9). In this low channel SNR regime, we have a closed form solution for the Cauchy 
case: 



UJr 



tan~ 



9 



(19) 



0" 2 7 " 

If we let 9 — > oo in (19), we get Uq = ir/9 which agrees with Theorem 2. 

3) Laplace Sensing Noise: In this case, we have (p n (u>) = 1/(1 + b 2 u 2 ) and b 2 := cr 2 /2. 
Substituting this in D(u) and letting (3 = u9, and differentiating D(j3/9) with respect to f3, 
equating to and simplifying we get, 

1 



G L (f3) 



b 2 



-i 2 



# 2 ^ 



tan 



a 







(20) 



with a := [1 + (cr^/P T )]- It can be easily verified that equation (20) has a unique solution in 
f3 E (0, 7r) as shown in Appendix 3. The /3£ that solves (20) can be found numerically and 



= 131/0. 



4) Uniform Sensing Noise: For the uniform sensing noise, we have (p n (u>) = sm(ua)/ua, 
where a 2 = a 2 /3. Substituting (p n (u>) in (11) and letting j3 = uoa for convenience we have 

[l-cos(f)] 



[a/3 2 csc 2 (/3) - 1] 
where C(f3) := 1 /[a/3 2 esc 2 (/3) - 1]. Writing D'(/3) = we get 



cos 



[a/3 2 esc 2 (0) - 1] - a/3 



2a 
"0 



tan 



9_ 

2a, 



6 



esc 2 (13) [1-/3 cot 







(21) 



(22) 



with a := [1 + (a 2 /P T )]. Theorem 1 does not apply for the uniform sensing noise. However if 
9/ a > 2, then using C(j3) > C((3+k7c), k = 1,2, ... , and using the periodicity of [1— cos((39 /a)], 
we can show that /3{j G (0,na/9}. Following similar arguments to the Laplacian noise case, it 
can be shown that there is only one stationary point in (0, ira/9] which corresponds to the global 
maximum. The /3{j that solves (22) can be found numerically and therefore, = /3^/ a. On the 
other hand if 9/a < 2, multiple local maxima are possible in (3 G (O,7ra/0] and (22) can have 
multiple solutions. In this case, that /3{j which yields the largest value for D(/3) in (21) should 
be chosen. 



November 16, 2010 



DRAFT 



11 



C. Per-sensor Power Constraint or high Channel SNR 

We now consider the DC under the per-sensor power constraint. In this setting, as L — > oo, 
Pr — > oo which makes (cr^/Pr) — > 0. Therefore the DC for the per-sensor constraint when L is 
large becomes: 

Equation (23) can also be interpreted as the DC when a 2 = for any finite L. In what follows, 
we characterize oo* in this per-sensor constraint regime, which effectively amounts to the removal 
of (a 2 /P T ) from (11). In this case there is not necessarily a uj* that attains the maximum in 
(12). Our first result reveals that (23) can be made large by choosing to sufficiently close to 
zero when rii are Gaussian, and yields an interesting relationship between the DC and the Fisher 
information. 

Theorem 3. When rii are Gaussian, 

supP pspc (w) = — = lim P pspc (w) (24) 

Proof: We begin with the inequality [1 — cos(u9)] < u 2 9 2 /2. Consider [20, eqn (1)], which 
using the fact that (f n (uj) is real-valued, reveals tp 2 n {u) < (1 + (p n (2u))/2. Using these two 
inequalities we can write the following: 

1 > L 1 :^- (25) 



D pspc (u) ~ 2u 2 V 2 n (u)9 2 ■ 
Now from [20, eqn (2)] with the fact that y? n (w) is real- valued, we have: 

[l^jgngaQj 1 
2u*<p*( U ) - J K } 

where J is the Fisher information of with respect to a location parameter [21, eqn (8)] (i.e., 

the Fisher information in xi about 9). Combining (25) and (26) we have: 

^M< J= i_ (27) 

9 2 a 2 n 

where the equality follows from the fact that for Gaussian random variables the Fisher infor- 
mation is given by the inverse of the variance. Now, we also see that using l'Hospital's rule on 
(23), lim^o Ppspc(w) = 9 2 /a 2 , which shows that the inequality in (27) can be made arbitrarily 
tight establishing sup w>0 D pspc (cu) = 9 2 /a 2 n . ■ 
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The proof of Theorem 3 also reveals an interesting inequality between the DC and the Fisher 
information, which of course is related to the Cramer-Rao bound for unbiased estimators. So 
for the per-sensor power constraint case with Gaussian noise, to should be chosen as small as 
possible for the best performance and it does not depend on the value of 9. 

For the Laplacian case, the solution is similar to the Gaussian case. It can be easily verified that, 
with (ct^/Pt) = 0, D' pspc (u) < over to e (0,n/9). This means that D pspc (u) is monotonically 
decreasing with to which implies that u should be chosen arbitrarily small. 

On the other hand, when n» are Cauchy distributed, then (p n (co) = e~ lu . Substituting in (23) 
and using 1' Hospital's rule we observe that lim w _>o -D pS p C (w) = for Cauchy sensing noise. This 
implies that, for the Cauchy sensing noise with per-sensor power constraint, smaller values of 
u should be avoided for reliable detection to be possible. 

D. Analysis of the DC for Non-homogeneous Sensors 

Consider now the case where rii are independent with non-identical distributions. This could 
occur if rii have the same type of distribution (e.g. Gaussian) with different variances. Letting 
(Pm(u>) = E[e jojni ], the DC in (10) becomes 



which is now a function of L unlike in (11), and reduces to (11) if op ni (co) = y?n(w), as in the 
iid case. We now study the conditions on the variances of := var(nj) for lirn^oo D L {uS) = 
for all cu > 0. When this asymptotic DC is zero for all cu > 0, the interpretation is that there is 
no suitable choice for cu > 0. The following result establishes that if the sensing noise variances 
are going to infinity, the asymptotic DC is zero for all to > 0, indicating a regime where reliable 
detection is not possible. 

Theorem 4. Let = ¥n(ciu) for some CF tp n (uj) where n has an absolutely continuous 

distribution. Suppose also that Hindoo <jj = oo. Then lim^oo Dl(lu) = for all u> > 0. 

Proof: Clearly the denominator of (28) is bounded between (<t^/Pt) and (1 + a^/Br). 
Therefore, it suffices to show that L^ 1 J2i=i VnAu) = L~ l J2i=i Vn{&iw) — > as L — > oo. Since 
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n has an absolutely continuous distribution, Hindoo <p n (x) = 0, and because lim^oo a { = oo, 
it follows that lim^oo <^ n (<7j0;) = for oo > 0. From [22, pp. 411] we know that if a sequence 
satisfies lim^oo aj = then its partial sums also satisfy lim^^ L _1 £ i=1 a$ = 0, which gives 
us the proof when applied to the sequence (p n (aiU>). ■ 
If, instead of of — > oo as « — > oo, the variances of are bounded, we can show the existence 
of an oj > for which lirn^oo D L (oo) > which is done next. 

Theorem 5. Let var(nj) exist for all i and o max := supj(var(n i )) 1 / 2 be finite. Then any < 
u < v^/omax satisfies Hindoo D L (cu) > 0. 

Proof: To show lim^oo D L (ui) > for w > 0, it suffices to show that L~ l J2i=i VnXu) > 
for uj > 0. From [23, pp. 89] we have tpr H {u) > 1 — crfu) 2 /2 for any CF with finite variance. 
Therefore, £ _1 £f =1 wM > 1 - £f=i W 2 > 1 ~ ^Lx^V 2 > where the last 
inequality holds provided that u < V2/a max . ■ 
This shows that if the noise variances are bounded, there exists (a small enough) u that yields 
a strictly positive asymptotic DC, establishing that there is a choice of uo that enables reliable 
detection. 



VI. Fading Channels 

Suppose that the channel connecting the i th sensor and the FC has a fading coefficient hi := 
Ihile^ 1 normalized to satisfy E[|/ij| 2 ] = I. If the sensors do not know or utilize their local 
channel information, and the fading has zero-mean (E[/ij] = 0), then the performance over 
fading channels is poor because the DC in (10) becomes zero due to law of large numbers and 
reliable detection is not possible. On the other hand, if the i th sensor corrects for the channel 
phase before transmission, using local channel phase information, the received signal under the 
TPC becomes 




Jo*. i Pt 



L 



Hl = ■ V — ^2\hi\e jwni +v , (29) 

=i 

where we focus on the iid sensing noise case to highlight the effect of fading even though the 
non-homogeneous case can also be easily pursued. The phase correction does not change the 
constant power nature of the transmission. By calculating the expectations in (10), for the signal 
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cos(u6)} 



(30) 



We see that in case of fading, the term y? 2 (w) is scaled by the factor (E[|/ij|]) 2 in the DC 
expression. Since E[|/ij| 2 ] = I, using Jensen's inequality, the factor (E[|/i;|]) 2 < 1 unless \hi\ is 
deterministic in which case it is one. Comparing (11) and (30) we have, with (E[|/ij|]) 2 < 1, the 
numerator of (30) is decreased and the denominator of (30) is increased, leading to a reduction 
in DC and thus fading has a detrimental effect on the detection performance, as expected. 

Note that if the optimization of the DC is desired in the fading case, the factor (E[|/i;|]) 2 in the 
denomenator of (30) affects the optimum u value. Theorem 1 can be proved for the fading case 
as well with C(u) := 2(E[|/i;|])V 2 (w)/[l - (E[|/i;|]) 2 y9 2 (w) + a 2 /P T ] which is still decreasing 
with cj if (p n (u>) is. Therefore the conclusion of Theorem 1, namely, cu* G (O,7r/0), does not 
change. The procedure to find the u* under the TPC for Gaussian, Cauchy and Laplacian is the 
same as described in Sections V-Bl, V-B2 and V-B3 respectively. The equations (17), (39), (40) 
and (42) remain valid with the exponentials in these equations scaled by the factor (E[|/ij|]) 2 . 
The equations (20), (44) and (45) for the Laplacian case also remain valid except the term 1/a 
in (20) scaled by (E[|/i;|]) 2 . 

We note that if sensors have imperfect knowledge of the phase, \hi\ will be replaced by \h i \e j<t ' i 
where 0; is the phase error. Clearly this error can also be subsumed in (29) as replacing uorii 
with bjrii + 4>i which changes the sensing noise by a term independent of u. This establishes 
the interesting fact that phase error over fading channels can be treated as a change in sensing 
noise distribution. 

VII. Asymptotic Performance and Optimization of u based on error exponent 

The error exponent in a distributed detection system is a measure of how fast the probability 
of error goes to zero as L — > oo. Mathematically error exponent is defined as: 



Large deviation theory [24], [25] provides a systematic procedure to calculate the error exponent 
which is briefly reviewed next. Let Y L be a sequence of random variables without any assumptions 



— lim 



log AM 

L 



(3D 
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on their dependency structure and let M(t) = lim i _ i , 00 (l/L) \ogE{e tYL } exist and is finite for 
all t e R. Define 

(32) 



e(z) = — lim — logPr[Fi < z\ , 



where z is the threshold and Y L is the test statistic of a detector. Gartner-Ellis Theorem [24, pp. 
14] states that e(z) in (32) can be calculated using, 



where 



e{z) = sup[tz - M(t)} , 
teR 



M(t) = lim -logE{e tyL }. 



(33) 



(34) 



We will now use the Gartner- Ellis Theorem with Y L replaced by Z L (u>) in (9) and z = 0. Letting 
M u (t) := lim L ^ 00 (l/L) logE{e* Zi ^)}, and e u (z) = sup teR [tz - M u (t)] we have the following 
theorem which relates the error exponent to the CF ip n (u)) of the sensing noise distribution. 



Theorem 6. For the detector in (5), the error exponent in (31) is e u (0) = — inf t6fi M u (t) where 

' 7 k=i 

+ 



M u {t) = log 



Iq ( 2 sin 



t) +2^4(2 



:su) | — \t I v9 n (fcw)cos I fc I - - — 

tV2(l - COSM)) 



2Pt 



(35) 



where I kit) is the modified Bessel function of the first kind. 



Proof: Please see Appendix 4. ■ 
It is well known that the function M u (t) is convex in t [24]. Therefore the supremum in (33) 

can be found efficiently for z = 0. The t* that maximizes (33) satisfies M u (t*) = which can 

be found by convex methods with geometric convergence [26]. 

In addition to the error exponent, it is also possible to approximate P c (u) using the function 

e u (z). In fact Bahadur and Rao [25, pp. 10] have proved that this probability can be approximated 

using the error exponent and is given by: 

1 



-Le u (0)(l+o(l)) 



(36) 
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as L — > oo and a 2 := [^l;(0)] 2 /[e2(0)]. The quantities 6^(0) and e'^(0) are the first and second 
derivatives of e u (z) at z = respectively, and can be calculated from the following equations 
[26, pp. 121]: 

e u (0)=t*, (37) 
4(0) = ^y- (38) 

The error exponent given in Theorem 6 is a function of lu and let us denote it by for 
convenience. It will be illustrated in Section IX that the values of u that minimizes P e (co) is 
closely predicted by the value obtained by maximizing D(u) or e w . We will also examine in the 
simulations in Section IX how accurately (36) can be used to approximate P c (oo) for finite L. 

VIII. Non-Gaussian Channel Noise 

We have so far assumed that the channel noise as Gaussian. However, we verified that the 
detector in (5) works well even if the channel noise is mixed Gaussian, uniform or Laplacian. 
The channel noise distribution will only affect the error exponent through the second term in 
(35). Using this, the effect of different channel noise distributions we considered are briefly 
sketched below. 

We considered the case of mixed Gaussian having two different variances drawn from a 
Bernoulli distribution. Let p be the probability that the samples drawn from the mixture have 
variance a 2 o and p x = l—p be the probability corresponding to cr^ and let a 2 vi > a 2 Q . In this case, 
we found that the error exponent is affected only by the larger variance in the mixture. While us- 
ing Gartner-Ellis Theorem to calculate M u (t), the second term in (35) for the mixed Gaussian be- 
comes lim L _1 log \p exp (t 2 a 2 (l - cos(lu6))/2P t ) + p 1 exp (t 2 aL (1 - cos(u;0))/2P T )l and 
this limit evaluates to \t 2 a 2 Vi {\ — cos(o;^))/2Pt] which proves that only the larger variance o 2 x 
in the mixture affects the error exponent. 

For the uniform channel case, interestingly we found that the second term in (35) evaluates to 
and thus proving that the error exponent is not impacted by the uniform channel noise. We do 
not include the straightforward derivation due to lack of space. We will discuss the performance 
of the mixed Gaussian and Laplacian cases in Section IX-F. 
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IX. Simulations 

We define the sensing and channel SNRs as p s := 2 /a 2 , p c := Pt/ct 2 and assume Pi = P = 
0.5 throughout. Note also that p = Pt/L is the power at each sensor as defined in Section II. 

A. Effect of ou on Performance 

We begin by comparing the optimized uo values using D(u), e u and P e {u) for the TPC. The 
values of cu* > obtained by maximizing the error exponent e w and the DC D(uj) were found 
to be very close over the entire range of P T . Figure 1 shows the plots of D(u),s w , and P e (ui) vs 
uj for Gaussian and Cauchy sensing noise distributions where the P e (u) plot is obtained using 
Monte-Carlo simulations. The different u* values in Figure 1 correspond to the best ui values 
obtained by optimizing D(u), e u and P e (u) respectively. It is interesting to see that the lu* 
that minimizes P c (u) is very close to that which maximizes D(uo) and e w . For Laplacian and 
Uniform sensing noises (not shown), the same trends were observed. 

Figure 2 shows the performance under per-sensor power constraint with large L. It is observed 
that smaller to yields better error probability. This agrees with our findings in Section V-C where 
it was shown that D pspc (lu) can be made larger by choosing u > arbitrarily small. Since 
both Figures 1 and 2 verify that the choice of u based on minimizing P e (u>) can be closely 
approximated by that which maximizes D(uo), in all subsequent simulations, we have used the 
cu* values obtained by maximizing D(oo). 

B. Comparison against MAF and MDF Schemes 

In Figure 3, the proposed scheme is compared under the TPC with the MAF and MDF schemes 
which have been shown in [8] to outperform conventional amplify-and-forward (AF) and detect- 
and-forward (DF) schemes. We observe that the proposed scheme outperforms MAF when p s > 
4 dB, and MDF for the entire range of p s . The same trend was observed when L is increased 
to 90 with an improvement in the detection error probability. The ML performance shown was 
obtained by the Monte-Carlo implementation of the ML detector and is computationally complex, 
but serves as a performance benchmark. Figure 4 shows the P c performance versus L under the 
TPC. Clearly the proposed scheme outperforms the AF, DF, MAF and MDF schemes consistently 
since p s = 15 dB. 
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The proposed scheme requires the fine tuning of the transmission phase parameter u> either 
through optimizing the deflection coefficient or the error exponent. However, it should be noted 
that similar type of fine tuning is also required in the competing schemes such as the MAF or 
the MDF. We note that the proposed scheme is inferior to MAF at low sensing SNRs (p s < 
4 dB). On the whole, the benefits of constant modulus signaling and improved performance at 
higher sensing SNRs make the proposed approach a viable alternative. 

C. Total Power Constraint: Different Noise Distributions 

For the Total Power Constraint, Figure 5 shows that Cauchy sensing noise results in better 
performance when p s is low, and worse when p s is high compared with other sensing noise 
distributions. This agrees with the fact that D(u*) is smaller for Cauchy sensing noise when p s 
is high than other distributions and vice versa when p s is low. When p s is moderately high, we 
observe that Gaussian, Laplacian and Uniform distributions have identical performance if p c is 
very low for a wide range of L as illustrated in the Figure 5. We found numerically that the 
similarity of the P e (u>) curves under different sensing noise distributions was also reflected in 
the corresponding D(oS) values where they were also verified to be similar. 

Figure 6 compares the performance of the proposed scheme in the presence of Rayleigh flat 
fading between the sensors and the FC against without fading with the Gaussian sensing noise. 
Clearly, fading has a detrimental effect on the detection performance as argued in Section VI. 
It is also observed that, in the presence of fading, P c is not as sensitive to the increase in p s as 
that of the no fading case. 

D. Error Exponent 

Figure 7 depicts the error exponent of the proposed scheme under the PSPC and illustrates 
its improvement with increase in p s for all the sensing noise distributions. Recall that al has no 
effect on the error exponent for the PSPC case since (o^/P T ) — > in (35). It is interesting to 
see that Cauchy sensing noise has a better error exponent than Gaussian, Laplacian and Uniform 
sensing noise distributions when p s < 4 dB while it is worse when p s > 4 dB. The error exponent 
with Gaussian sensing noise is better than that of Laplacian noise when when p s > 7.5 dB and 
the uniform distribution has a better error exponent than other sensing noise distributions when 
p s > 4 dB. The error exponent of the proposed scheme is compared with the error exponents of 
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MAF and MDF schemes which were only derived for the Gaussian case (please see equations 
(24) and (25) in [8]). It is seen that, for the PSPC case, the MAF scheme (whose error exponent 
is p s /S) and the proposed scheme with optimum u have identical error exponents leading us to 
conjecture that sup^[— m.i teR M u (t)] = p s /8 when rij are Gaussian. The MDF error exponent is 
inferior compared to MAF and the proposed scheme. 

Figure 8 shows the error exponent under the TPC with p s = dB. In this scenario, Cauchy 
sensing noise has the best error exponent since p s is low. This concurs with the fact illustrated 
in Section IX-C that the DC of Cauchy is better at lower values of p s than other distributions 
and this was justified by the simulation results as shown in Figure 5. We found that when 
p s is increased, Cauchy becomes inferior to other noise distributions. For all the distributions, 
increasing p c results in an increase in the error exponent which becomes a constant beyond 
p c = 15 dB. This is because, for a given p s , increasing p c combats the effect of channel noise, 
thereby improving the error exponent. However, the effect of sensing noise can not be overcome 
by increasing p c indefinitely. This can be seen from (35) as well where the second term vanishes 
while the first term remains even for large P T . For the Gaussian case, we derived the error 
exponent of the MAF scheme under the TPC as £ MAF = 2 /8[a 2 n + {o 2 v {a 2 n + P Pi6 2 )/P T )]. If 
Pt — > oo, this reduces to p s /8 for the PSPC case. It is seen that with p s = dB, the MAF 
scheme is better than the proposed method when p c < 15 dB. However, under the TPC, the 
error exponent of the proposed scheme was found to beat the MAF scheme when p s > 4.5 dB 
and an example plot is shown in Figure 8 for p s = 10 dB. This crossover between the MAF and 
the proposed schemes is also reflected in their respective P c performance curves approximately 
around the same p s value (please see Figure 3). However, if p c is increased beyond 15 dB, we 
see that the error exponents of both the schemes become very close. 

E. Approximations of P e (ou) through e u (z) 

Equation (36) provides an approximation of P c (u) based on the error exponent. The expression 
in (36) is found to match with the simulations when p c > dB and p s > -5 dB. Figures 9 and 
10 elucidate this behavior for Gaussian sensing noise distribution. Similar trends were observed 
for the other sensing noise distributions as well but are not shown due to space constraints. 
When L is small, the gap between theory and simulation is significant as shown in Figure 9. 
This can be explained by the o(l) term in (36). Accordingly, when L is increased to about 40, 
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we see the theory and simulation curves merging as shown in Figure 9. Figure 10 shows that 
when p s is moderately high, smaller L is required to get the performance match between theory 
and simulation. 

From the various simulation plots in Figures 1, 5, 7, and 8, we see that the proposed scheme 
is robust in the sense that it works very well for a variety of sensing noise distributions including 
the impulsive Laplacian distribution and the Cauchy distribution which has no finite moments. 

F. Non-Gaussian Channel Noise 

Figure 11 shows the error exponent plot for the case where a 2 = 0.25, po = 0.80, = 
4, p! = 0.20 (note that the effective channel noise variance is: a\ = Po<t% + Pi&vi = !)• We 
see that the error exponent of mixed Gaussian with a% eS = 1 is worse compared to that of the 
Gaussian with al — 1 case. This is because, in the mixed Gaussian case, the error exponent is 
a function of the larger variance of cr^ = 4. 

Figure 12 shows the performance of the proposed detector with Laplacian channel noise against 
the Gaussian channel noise when the sensing noise is Gaussian. We note that when sensing SNR 
p s is moderately high, the impulsive Laplacian channel noise is worse compared to Gaussian 
channel noise. 

X. Conclusions 

A distributed detection scheme relying on constant modulus transmissions from the sensors is 
proposed over a Gaussian multiple access channel. The instantaneous transmit power does not 
depend on the random sensing noise, which is a desirable feature for low-power sensors with 
limited peak power capabilities. The DC of the proposed scheme is shown to depend on the 
characteristic function of the sensing noise and optimized with respect to oo for various sensing 
noise distributions. In addition to the desirable constant-power feature, the proposed detector is 
robust to impulsive noise, and performs well even when the moments of the sensing noise do 
not exist as in the case of the Cauchy distribution. Extensions to non-homogeneous sensors with 
non-identically distributed noise are also considered. It is shown that over Gaussian multiple 
access channels, the proposed detector outperforms AF, DF and MDF schemes consistently, and 
the MAF scheme when the sensing SNR is greater than 4 dB. The proposed detector is shown 
to work with the non-Gaussian channel noises as well. The error exponent is also derived for the 
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proposed scheme and large deviation theory is used to approximate P e (uj) for large L. It is shown 
that while the DC has a simpler expression for the purpose of optimizing ui, the probability of 
error approximation based on (33) is shown to be an accurate indicator of detection performance 
for all distributions and moderate number of sensors. The effect of fading is also considered, 
and shown to be detrimental to the detection performance. 

Appendix 1 : Gaussian Sensing Noise 

First we note that G G (0) = (a — 1) > since u 2 > and G G (ir) = — oo. Since Gq(P) is 
continuous, (17) has at least one solution. To show that this solution is unique, consider the first 
derivative: 



2a I — sec 2 ( — 



I an I — 



(39) 



2 V 2 . 

Now, using tan(/3/2) > (3/2 and sec 2 (/3/2) > 1 + (/3 2 /4) for p G (0,tt), we get the following 
upper bound: 



6 2 



2 U ^ -ap( 1 + ^- ] -n I 



(40) 



Since a 2 > we have a > 1. Recall that j3 e (0, n), and the rhs of (40) is always negative. It 
follows that Gq(P) is monotonically decreasing over j3 e (0,7r) and (17) has a unique solution 
which corresponds to the global maximum of D(/3/9). 

Appendix 2 : Cauchy Sensing Noise 



The first derivative of D((5/9) is given by, 



D 



sin 



(cte e P — If 



a — e » M —a tan I — 

p V 2 



Since the first term on the rhs of (41) is non-zero for (5 G (0, it), we need to solve 



Go(/3) :=a-e-^-^atan (| 



0. 



(41) 



(42) 



First we see that G c (0) = (a — 1) > and G c (vr) = — oo which implies that there is at least 
one solution to (42) in j3 e (0, 7r) as G c ((5) is continuous. The second derivative of Gc(P) is 
given by 



Gc(/9) = 



1 e 2 



7 « 2 m. //? 

H — — sec — tan — 
\2 \2 



(43) 



Clearly, G c {0) < for (3 G (0, tt) which establishes that G c (/?) is concave. Therefore, (42) has 
a unique solution which corresponds to the global maximum of D(/3/9). 
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Appendix 3 : Laplacian Sensing Noise 

First we note that G L (0) = (1 - (1/a)) > if a 2 > and G L (n) = -oo. This means that 
(20) has at least one solution. The first derivative of C7 L (/3) is given by, 
2b 2 



2/3 (1 + ^ 2 



P + P 3 ) -c 2 (|) + 2(l + 3^ 2 )tan(§ 



• (44) 



Now, using tan(/3/2) > (3/2 and sec 2 (/3/2) > l + (/3 2 /4) over /3 G (0, tt) in (44) and simplifying, 
we get the following upper bound: 

b 2 



G'M<- W [(6 2 + 8b 2 )f3 3 + b 2 f3 5 } 



(45) 



Clearly, for j3 e (0, 7r), the rhs of (45) is always negative which implies G' L (/3) < 0. It follows 
that is monotonically decreasing over j3 e (0, 7r) and (45) has a unique solution which 

corresponds to the global maximum of D(/3/9). 

Appendix 4: Proof of Theorem 6 

We use the Gartner-Ellis theorem from large deviation theory [24, pp. 14] to calculate the 
error exponent. To this end, we need to calculate M^it) in (34) and substitute into (33). 

M w (f) = lim ^logE{exp[tZ L ]} 



= lim — log E <( exp 



logE | exp 



[J^ (uj6\ ( uj9 tt\ fT 
t n 2 sin ^— J cos [un t -- + -j +] J— z 

t 2 a 2 (l -cos(u9)) 



. . u>6\ ( oj6 tt 
2t sin | — I cos \urii- — + - 



2P T 



(46) 



From [27, pp. 376], we have the Fourier series expansion of the periodic function e pcos(M ) as, 

oo 

e pcos(u) = I (p) + 2 Up) cos(ku) (47) 



k=l 



Using the equation (47) in (46) with p = 2tsm(u6/2) and u = (urn* — uo9/2 + n/2) and then 
applying the expectation on the resulting summation, we get M u (t) as in (35). 
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Fig. 1: Total Power Constraint, D(u),e u , P e (u) vs u: p s =10, 15 dB, p c =-10 dB, L=20 
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Fig. 2: Per-sensor Power Constraint, Gaussian, P e (u) vs u: p s =-10 dB, p=10 dB 
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Fig. 4: Total Power Constraint, P c vs L: p s =l5 dB, p c =0 dB 
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Fig. 5: Total Power Constraint, P c vs L: p s =5, 15 dB, p c =-10 dB 
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Fig. 9: Gaussian Sensing Noise: P e vs p c : p s =0 dB, L=10, 20, 40, 60 
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Fig. 10: Gaussian Sensing Noise: P c vs p c : p s =l0 dB, L=5, 7 
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-9 — Mixed Gaussian 



-a — Gaussian, o =1 
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Fig. 11: Mixed Gaussian channel noise: e u vs p c : o^ o =0.25, p = 0.8, 0^=4, Pi = 0.2, = 1, 
p s =0 dB. 
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Fig. 12: Gaussian Sensing Noise: P c vs p s : Px=- 12.22 dB, L=60 
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